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Abstract 



It is known that there is a transition from photon antibunching to bunching in 
the resonance fluorescence of a driven system of two two-level atoms with dipole-dipole 
interaction when the atomic distance decreases and the other parameters are kept fixed. 
We give a simple explanation for the underlying mechanism which in principle can also 
. be applied to other systems. 

m ■ PACS numbers 42.50.Ar, 42.50Fx 

o ■ 

00 ■ 1 Introduction 

^ . 

, Bunching means that photons emitted by a driven system in steady state have a 

tendency to arrive in pairs or larger groups at a detector rather than uniformly 
^ ■ distributed in time. More precisely, right after a photon emission the probability 

^ . density for emitting another photon is larger than for a uniform distribution of 

Q^! corresponding emission rate. For antibunching this probability density is smaller 

than for a uniform distribution, and it means that the photons repel each other. 
K/i , For a driven two- level system it is known that one has antibunching [|l], |^, |^. 

' This is intuitively very simple to understand since after an emission the system 

is in its ground state and it requires some time to acquire enough population 
of the excited state for a next emission. For two independent, noninteracting, 
two-level atoms one also has antibunching, although not quite so pronounced as 
for a single atom. This was investigated experimentally in Ref. Q. 

For a system of two two-level atoms with dipole-dipole interaction it is known 
from studies of the master or optical Bloch equations for the system that the prop- 
erties of the resonance fluorescence may change considerably when compared to 



that for a single system p[-|26]. In particular there is a transition from anti- 



bunching to bunching when the atomic distance becomes small and the other 



parameters are kept fixed (see e.g. Ref. [[T^]). The state space of the system is 
four- dimensional and the density matrices have 16 components so that the corre- 
sponding Bloch equations require diagonalization of a 16 x 16 matrix. This makes 
an intuitive understanding of this transition from antibunching to bunching not 
obvious. 
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It is the aim of this paper to elucidate the underlying reason and give a 
simple explanation for the appearance of bunching for small distances in a driven 
system of two two-level atoms with dipole-dipole interaction. We will trace the 
phenomenon to two causes. One is the form of the steady-state density matrix 
of the system. The other is the density matrix of the system right after an 
emission, called the reset matrix which depends only on the state prior to 
the emission. In it the steady-state ground-state population has disappeared 
since it does not contribute to an emission, while the steady-state populations of 
the higher levels have been transferred to levels one step lower, in proportion to 
their respective decay constants, and then normalized to 1. The new populations 
of the intermediate states then determine the probability density for the next 
emission. 

The transition from antibunching to bunching can be directly read off from 
the steady-state density matrix of the system without a lengthy calculation. For 
small atomic distances and for the Rabi frequency of the order of the Einstein 
coefficient A or a few times larger the ground-state population is very large while 
the excited states have small but similar populations, and thus there is only a 
small steady-state emission rate. The large ground state population does not 
contribute to a photon emission and therefore disappears after an emission, while 
the populations of the intermediate states move to the ground state and that of 
the highest state to the intermediate ones, in proportion to their respective decay 
constants and with ensuing normalization. Hence after a photon emission the 
ground and excited states have suddenly acquired populations of the same order 
of magnitude and thus a higher emission probability density than before. This 
means that the photons tend to come in pairs and that there is bunching. For 
weaker driving and small distances the same mechanism holds. In this case the 
population of the ground state is even larger and that of the highest excited state 
smaller than for the other excited states. But its population is still large enough 
so that after an emission, when the populations have moved one step down, the 
new population of the intermediate states is larger than before so that there is 
an increased emission probability and thus bunching. These considerations can 
be applied to quite general systems. 

Here we consider a system of two two- level atoms at a fixed distance r, in- 
teracting with the quantized radiation field and a classical laser field. Through 
photon exchange the radiation field mediates the r dependent dipole-dipole in- 
teraction of the atoms. The dipole and rotating-wave approximation is used 
throughout. Retardation effects are included. 

In Section 2 we briefiy review the photon-counting correlation function g{T). 
If g{0) > 1 one has bunching, if g{0) < 1 one has antibunching. In Section 3 we 
apply this to two two- level atoms and give an explicit expression for g{0) as a 
function of the atomic distance and the driving field. Bunching is explicitly seen 
for atomic distances about a quarter of a wavelength or less. 

In Section 4 we discuss the results and the simple mechanism responsible for 
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bunching for small atomic distances. It is also made clear why an analogous 
argument yields antibunching for two independent non-interacting atoms and for 
atoms sufficiently far away. 

For simplicity we consider, in the main part of the paper, coinciding atomic 
dipole moments and the same laser phase for each atom. In the Appendix the gen- 
eral case is considered, the corresponding Hamiltonian spelled out and the reset 
matrix given. We also outline the connection with the quantum jump approach 
[E^, p9|, p^, pOi which is equivalent to the Monte Carlo wave-function approach 



3ll] and to quantum trajectories For a recent review of this approach see 
Ref. [^]. For a system of two two-level atoms with dipole-dipole interaction we 



carry this approach over in a form convenient for simulations 



2 Photon-counting correlation functions 



We briefly review some well-known facts. As pointed out in Ref. ||3J] there is 
a minor difference between correlation functions which are based on the electric 
field operator and those which are based on the photon number operator. Here 
we employ the latter type correlation functions of second order. For simplicity we 
consider broad-band photon detections over all space. It is useful to distinguish 
clearly between correlation functions for ensembles and for a single trajectory. 
The latter involves a time rather than an ensemble average. 

Ensemble. Consider an ensemble of laser driven atomic systems in the state p 
at t = and denote by G{t2, ti; p)dtidt2 the relative number of systems for which 
in addition to a photon in [ti, ti+dti) also a photon in (^2, t2 + dt2) is detected. If, 
for a particular trajectory, we denote the number of photons detected in (t, t+ At) 
by A^*''^j (t, t + At) - for small At this number is either or 1 - then 

G{t2M\ P)dtldt2 = {N''-''^{ti,ti + dt,)N''^^{t2,t2 + dt2))ens ■ (1) 

Let us consider the sub-ensemble of systems which had an emission at ti and let 
us denote its normalized density matrix right after the emission by TZ{p{ti)). This 
we call the normalized reset matrix and it will be given explicitly for a two-atom 
system in the next section. We denote by I(t; p) the probability density for the 
emission of a photon at time t (not necessarily the first photon after t = 0) for 
initial density matrix p. With this one has 

G(t2, ti; p)dtidt2 = /(ti; p)dti /(ta - h; 'R{p{ti)))dt2 . (2) 

Letting ti ^ oo and keeping t = t2 — ti fixed the first factor on the r.h.s. 
goes to Jss and p(ti) to p^^, the steady-state emission rate and density matrix, 
respectively. Hence 

G{t) = hm G{t + r, t- p) = hJ{T- nip^) . (3) 
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The photon correlation function g{T) is now defined as 

«W.^^^<4^^. (4) 

-'ss-'ss -'ss 

It compares, for the steady state, the probabihty density for emission of a photon 
at a time interval r after a preceding emission with that of a uniform distribution 
of emission rate Igs- 

Single trajectory. We now consider a single system with its trajectory of 
photon emissions and define A^*''^j as before. At instances t'^ = mAt' until 
time T = MAt' one measures whether or not a photon has been emitted in 
it'^,t'^_^_j). Then the relative frequency of cases in which both in {t'^,t'^ + Ati) 
and {t'^ + T,t'^ + T + At2) a photon has been found is given in the limit At' — ^ 
and T ^ oo, using 1/M = At'/T, by 

At' 

G*-j(r)AtiAt2 = lim hm — J]iV*-j(C4 + Ati)iV*-j(4 + r,4 + r + At2) 

m 

= \im ^ [ &t' N^'^'^t'^t' + Ati)N''^\t' + T,t' + T + At2) . (5) 

By ergodicity this should be the same for each trajectory, and therefore one can 
take the ensemble average of the r.h.s. without changing anything. Using Eqs. 
(|lD and (I) one then obtains 

Gf-^j(r) = G(r) (6) 

so that both correlation functions coincide and similarly g^^^^r) = g{T). 

We also point out the well-known fact that if one observes photons with a de- 
tector of efficiency t] less than 1 then in Eq. both numerator and denominator 
are multiplied by rj and hence gi^r) is not affected by the detector efficiency. 

One has bunching if the relative number of cases, in which shortly after emis- 
sion of a photon a further photon is emitted, exceeds those for a uniform dis- 
tribution of frequency /ss- Thus bunching means ^'(0) > 1. Similarly one has 
antibunching if this number is less than for a uniform distribution, i.e. if (7(0) < 1. 



3 Bunching for two atoms 

We now turn to two two-level atoms with dipole-dipole interaction, driven by 
a laser tuned to the atomic transition frequency ujq. The corresponding Hamil- 
tonian is given in the Appendix. For simplicity we consider coinciding atomic 
dipole moments forming an angle with the line connecting the atoms and laser 
radiation normal to this line so that the laser is in phase for both atoms. The 
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Rabi frequency of the laser, denoted by f2, is then the same for both atoms. One 
can take Q to be real and positive. The general case is indicated in the Appendix. 

It is convenient to use the Dicke states [^] \g) = |1)|1), |e) = |2)|2), and 
|s) and \a) the symmetric and antisymmetric combinations of |1)|2) and |2)|1). 
These states play the role of dressed states for the atoms (cf. e.g. Ref. pl|), 




hlmC 



Figure 1: Dicke states and decay rates 



with decay constants A ± Re C (see Fig.|l|) where C is an r dependent complex 
coupling constant. It is given for the general case in Eq. (|TB|) of the Appendix. 
From Fig. |] it is seen that C —>■ for r oo, |Im| C ^ oo for r ^ 0, while Re C 
changes little with r. Retardation effects are included in the sense that C goes 
to its value for a static dipole-dipole interaction when c 
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iJ = 
= 7r/8 
= 7r/4 
i5 = 7r/2 
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Figure 2: Dependence of C on r 

The steady-state density matrix p^^ can be found from the Bloch equations 
and is known in the literature, see e.g. Ref. |]12|. One can also directly employ 



5 



the Bloch equations in Eq. (^) of the Appendix and put p = p'^^ = 0. In the 
Dicke basis one obtains for the diagonal elements 

,3 _ (A^ + n^f + A^ReC i2A + Re C) + (Im C)^ 
Pgg ~ ]Y 

ss V I " " y ss ss 2_ ('7\ 

rss jy- 5 raa ree V / 

with the normalization factor 

N = {A^ + 2Q^y + A^ReC{2A + ReC)+A^{lmC)\ (8) 

We also need the diagonal elements of normalized reset matrix, the density 
matrix right after an emission. Due to an emission the populations of the excited 
states in the Dicke basis move down one step to lower levels in proportion to 
their decay constants and the previous ground-state population disappears since 
it does not contribute to an emission. Normalization is then achieved by dividing 
by the trace, tr(.). This gives 

{g\n{pn\g) = {{A + ReC)pZ + {A-ReC)pZ} /tT{.) 
{s\n{pn\s) = {A + ReC)pZ/ti{.) 

{a\n{pn\a) = {A-ReC)pZ/tT{.) (9) 

and (e|'^(p****)|e) = 0. The complete reset matrix is given in the Appendix. 

One can immediately draw the following conclusions from these expressions. 

(i) For small atomic distance, kor < 2, ImC and N become very large. Hence, 
both for weak and stronger driving, the steady-state population p^^ of the ground 
state becomes much larger than that of the excited states, and thus the steady- 
state emission probability is small in this case. 

(ii) For strong driving, ~ A, the ratios of steady-state populations of the three 
excited states are of equal order of magnitude for all atomic distances (since Re C 
does not vary much and drops out). 

(iii) Right after an emission, the (large) steady-state ground-state population 
is discarded, the populations of \s) and \a) are transferred to \g) in the reset 
matrix and that of |e) to \s) and |a), all in proportion to their appropriate decay 
constants. Hence, after an emission and for f2 ~ A, the ground-state population 
and that of the first excited states have become of similar magnitude (see Fig. || 
for a qualitative description). 

(iv) After an emission therefore, for small atomic distance and for ~ A, the 
population of the two first excited states has increased in relation to the ground- 
state population. Therefore the probability density for the next photon right after 
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Figure 3: Relative populations before and after an emission. 



an emission is higher than the steady-state emission rate. This means bunching. 

This argument for bunching can be extended to weak driving and small dis- 
tances as follows. From Eq. (J^) one sees that in the steady state only the 
numerator of the ground-state population contains ImC. Since the latter in- 
creases rapidly for decreasing r, as seen from Fig. ^ the ratios of the populations 
of the excited states with that of the ground state approach 0, while the ratios 
among the excited states do not change. After a photon emission the upper 
populations move downwards and the previous large ground-state population is 
discarded. Hence again, after a photon emission the ratios of the populations of 
excited states and ground state have increased compared to the steady state if the 
atomic distance is sufficiently small, and this means a higher emission probability 
density, i.e. bunching. 

These observations will now be made quantitative. Since Iss is obtained from 
the level population multiplied by their decay constants A±ReC, one has 

/ss = {A + ReC)pZ + {A-ReC)p::, + 2ApZ . (10) 

Hence the normalization constant tr(.) in Eq. (H) is Iss- For small atomic distance 
Iss becomes very small, due to the small population of the excited states. This 
can be attributed to the detuning due to the level shift hlmC (see Fig. jT]). 

For g{0) in Eq. (^) one needs /(O; 7?.(p*^'^)), the probability density for a new 
emission right after an emission. This is obtained in a similar way as Igg, 

l{0;n{p'')) = {A + ReC){s\n{p'')\s) + {A - ReC){a\n{p'')\a) + 2A{e\n{p'')\e) 

= 2{A' + {ReCY}pfJIss (11) 

by Eqs. and (p!0|). One could also have used Eq. (|T9|) of the Appendix. From 
the behavior of Re C it follows that /(O; TZ{p^^)) is of the same order of magnitude 
for all atomic distances. This fact is immediately understood by the observations 
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(ii) and (iii) above. From Eqs. (H) and (0) - (|TT]) one finally obtains 



9(0) 



/(0;7^(p-)) 



2^2 



1 + 



{2Q^ + A'^ + AReCf 



(12) 



Since Iss becomes small for small r while /(O; TZ{p^^) does not change much with r 
one has 5^(0) > 1 for small atomic distances. In the last expression the first factor 
approaches 1 for small atomic distance since ReC goes to A, while the second 
factor grows with ImC. In particular, for weak driving the terms involving Q 
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Figure 4: 
n = 0.9A 



Photon correlation g{0) as a function of r. (a) = 0.1 A, (b) 



can be neglected and one can read off Fig. |^ that one has bunching below an 
atomic distance of about a quarter of the optical wavelength. For strong driving, 
fl ^ A, bunching sets in when the atoms are slightly closer. For large atomic 
distance g{0) approaches 1/2 since C approaches 0. This recovers the result for 
two independent atoms. This is plotted in Fig. |. 



4 Discussion 

We have investigated bunching and antibunching in the resonance fluorescence of 
two atoms as a function of their distance and with their dipole-dipole interaction 
taken into account. Each atom was treated as a two-level system and the position 
of the atoms was kept fixed. The two-atom system was irradiated by a laser tuned 
to the transition frequency of the individual atoms. Retardation effects have been 
included. 

For a single two-level atom antibunching in the resonance fluorescence is well- 
known and well understood After emission of a photon the atom is in its 
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ground state and has to be pumped to the excited state before emitting a new 
photon. Hence the probabihty density for finding a second photon right after an 
emission is zero. 

For two independent, noninteracting two level atoms one of the atoms is in 
its ground state right after an emission while the other is unchanged. Therefore 
the probability density for a second photon right after an emission is half of that 
in the steady state. This means g{0) = 1/2. 

For two interacting two-level atoms the emission statistics depends on the dis- 
tance. For atoms far apart the interaction is negligible and one has antibunching 
as for two independent atoms. For small atomic distances one has bunching. The 
main purpose of this paper was to get a better understanding of this phenomenon. 

Let \g) = |1)|1) and |e) = |2)|2) denote the states where both atoms are in the 
ground and excited state, respectively, and let \s) and \a) be the symmetric and 
antisymmetric combinations of |1)|2) and |2)|1) (Dicke states). We first discuss 
the case of strong driving. Then, for small distances, the steady-state ground- 
state population is much larger than those in \s), \a) and |e) while the populations 
of the latter are of similar (though small) magnitude, as indicated in Fig. |^. The 
reason for the small population of the latter is easily understood through the level 
shift of |s) and \a) due to the dipole force (see Fig. |l]). The reason for the similar 
magnitude of the population of \s) and \a) with |e) has been attributed to two- 
photon processes connecting l^f) with |e) [^. Now, once a system has emitted 
a photon the population of |e) is transferred to |s) and \a) and the population 
of the two latter to \g), in proportion to the respective decay constants and with 
ensuing normalization. The previous population of \g) has disappeared since it 
does not contribute to the emission. Thus right after an emission the populations 
of I a) and \g) are suddenly of similar magnitude while before the emission the 
population of the ground state was much larger. Hence right after an emission 
the probability density for finding another photon has increased when compared 
to that preceding the emission, i.e. compared to the steady state emission rate. 
This means bunching for small distances and strong driving. 

For weak driving the mechanism is in principle the same. Although in the 
steady state the populations of the excited states now are no longer of similar 
magnitude the ground-state population increases with decreasing distance much 
faster than the population difference between the excited states. This means that 
the population of |e) is not too small when compared to the population of \s) and 
I a). Therefore after an emission, when the populations have moved down one 
step, the combined population of \s) and \a) has increased. This means a higher 
emission probability density than in the steady state, i.e. bunching. 

We have shown that when decreasing the atomic distance the transition from 
antibunching to bunching sets in at a distance of about a quarter of the optical 
wavelength, for weak driving slightly sooner than for strong driving. 

It is instructive to see why the same argument gives antibunching for two 
independent, non-interacting atoms. First, for strong driving, the two levels of 
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each individual atom are populated by approximately 1/2 so that the population 
of \g), \s), \a), and |e) are 1/4 each. Then, after an emission, the ratios of the 
populations of \g), \s), \a) and |e) are | : ^ : | : 0, as inherited from those of 

I a), |e), prior to the emission and in proportion to the decay rates. Thus the 
probability density for a next emission is only one half of that in the steady state. 
On the other hand, for weak driving the ground-state population is much larger 
than that of the excited states. Is this situation not similar to that of interacting 
atoms? Not quite, since although the populations of \s) and \a) are small and 
of similar magnitude, that of |e) is of the order of the product of the latter and 
therefore an order of magnitude smaller. Thus after an emission the ground state 
population is still much larger than that of the excited states, and there is no 
increase in the emission probability. 

The analysis can be carried over to the more general case where the dipole 
moments are not parallel and where the laser is detuned and its phase is different 
for the two atoms. The necessary tools are given in the Appendix. Also the case 
of degenerate upper level can be treated. The results are similar to those 
obtained above. 

To conclude, we have traced the appearance of bunching in the resonance fluo- 
rescence of a driven system of two two-level atoms with dipole-dipole interaction 
and at small distances to two causes, one the level populations of the steady- 
state density matrix, the other the change in the state right after the emission 
of a photon. A similar analysis can in principle also be applied to other systems, 
e.g. to a single atom in a three-level cascade configuration. 



Appendix 

We consider two atoms fixed at positions and each with two levels, |l)i and 
12)4,^ = 1,2, with energy difference huo. We define operators 5*^^ in the two- 
atom Hilbert space by = |2)jj(l| and = |l)jj(2|. The dipole moment 
of the 2-th atom is D^2 = j(l|X|2)j. For the laser we take zero detuning and 
EL(r,t) = Re{Eo e^C^^-'-'^o*) }. Making the usual rotatin g-wave approximation 
and going over to the interaction picture the interaction Hamiltonian becomes 

2 

Hi = Y,Y.^ [gSs^Ks e^("«-"^)* e''^-'-'5+ + h.c. ] + Hl , (13) 

1=1 k,s 

with the coupling constants 

laser part Hl = | Yl'i=i {^i^t + ^i^i} ^'^'^ Rabi frequencies i7j = | ■ 
Eq e^^^''^\ The operator Hj contains the dipole-dipole interaction of the two 
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atoms as seen from the Bloch equations or from the conditional Hamiltonian 
between emissions, as explained further below. In the above Power-Zienau for- 
mulation this interaction is due to photon exchange 

Reset matrix. The reset operation gives the state or density matrix right 
after a photon detection. In a basis in which the atomic damping is diagonal, 
as for the Dicke states, the diagonal states immediately can be written down, as 
in Eq. (P). For a general A^-level system the reset matrix has been derived in 
Refs. p7| , |30| . For a system consisting of two or more atoms the derivation has 
to slightly modified since in this case the field operator E appear with different 
position arguments. 

Let at time t the state of the combined system, atoms plus quantized radiation 
field, be given by |Oph) p (Oph|, i.e. the atomic system is described by the density 
matrix p and there are no photons (recall that the laser field is treated classically) . 
If at time t + At a. photon is detected (but not absorbed) the combined system 
is in the state 



FyUiit + At, t)|Oph) p (Oph|f/J (t + At, t)IP> 



(15) 



where i?> = 1 — |0ph)l/i(0ph| is the projector onto the one or more photon space 
(since At is small one could directly take the projector onto the one-photon 
space). The probability for this event is the trace over Eq. (0). For the state of 
the atomic system it is irrelevant whether the detected photon is absorbed or not 
(intuitively the photon travels away and does no longer interact with the atomic 
system). Hence after a photon detection at time t + At the non- normalized state 
of the atomic system alone, denoted by lZ{p)At, is given by a partial trace over 
the photon space. 



n{p)At = trph (lP>Ui{t + At, t)\%h) P {OphluUt + At, t)iP> 



(16) 



We call TZip) the non-normalized reset state |2^. Proceeding as in Refs. ||2^, ^ 

C21) ^2 pSi 



and using perturbation theory one obtains p6 

n{p) : 



I {CI, + C21) S^pSt + \ (C12 



+A {S^pSt + S^pSt) 
with the r dependent constants 

1 



(17) 



^12 1 ^12 



^) '-'12 



+ 



{kory i{kory 



^12 ' ^12 



3D 



12 > ■ 
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Here " denotes vectors normalized to 1, A is the Einstein coefficient, and r = 
r2 — ri. In the case of equal dipole moments one has C12 = C21 = C which was 
depicted in Fig. ^, with d defined by cos^ § = |(Di2, r)|^ /r'^D'12- The normalized 
reset state is 7^(p) = lZ{p) /tilZi^p). 

By Eq. (|T5|) the normalization of 7^(p) is such that trA7?.(p)At is the proba- 
bility for a photon detection at time t + At if the (normalized) state of the atomic 
system at time t is p. Hence one has for the probability density / of Section 2 

/(^;p(0)) = tr7^(p(^)) . (19) 

The laser field does not appear in the reset state, just as in the case of a single 



atom |3^ , since its effect during the short time At is negligible. By a simple 



calculation one checks that Eq. (^) can be written as 



7^(p) 



[a+]^ |Ci2 + C2M^ R+pRi \ \Ci2 + C2^il) R-pRl (20) 



where i?± = [S^ ± e'^S^) /V2 and if is the argument of C12 + C{2- From Eq. 
(|I8D one can check that A > \\C12 + Cgil- If p is a pure state, p = say, 
then i?±pi?^ are also pure states. This decomposition of 7^(p) is advantageous 
for simulations of trajectories. 

Conditional Hamiltonian and waiting times. In the quantum jump approach 



2^, |2^, 0, the time development of an atomic system is described by a 
conditional non-hermitian Hamiltonian ifcond, which gives the time development 
between photon emissions, and by a reset operation which gives the state or 
density matrix right after an emission. For a general A^-level system these have 
been derived in Refs. p7| , |5D| . The derivation of the former is adapted here to a 
system of two atoms. 



As explained in Refs. p8|, |30|, H^ond is of the general form -f^cond = 
Ha + Hl + r where F is an atomic damping operator. In a basis in which F 
is diagonal the diagonal terms are just the decay constants of the corresponding 
states. If these (dressed) states are known ifcond can immediately be written 
down. In this way one can obtain i^cond for parallel dipole moments in the Dicke 
basis. In the general case it is obtained (in the interaction picture) from the 
short-time development under the condition of no emission, i.e. from the relation 

1 - ^ i^condAt = {%h\Ul{MM%h) 

where the r.h.s. is evaluated in second order perturbation theory for At interme- 
diate between inverse optical frequencies and atomic decay times. In a similar 



way as for a single atom p8l BQ] one obtains for two two-level atoms [36 



i^cond = I [StS^ + 8^82) + C12 8^82 + C21 8^8^] + Hl (21) 
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with the r dependent constant Cij given by Eq. (0). Between emissions the time 
development is given by Ucondit, 0) = exp {—iHcondt/h} which is non-unitary since 
-f^cond is non-hermitian. The corresponding decrease in the norm of a vector is 
connected to the waiting time for emission of a (next) photon. If at t = the 



initial atomic state is then the probability Po{t) to observe no photon until 
time t by a broadband counter (over all space) is given by ||28|, ^7\, Q 



Po(t;|V'))=l|f/cond(t,0)|^)f , (22) 

and the probability density wi of finding the first photon at time t is 

w,{t;m = -j^Po{t;m. (23) 

For an initial density matrix instead of the expressions are analogous, with a 
trace instead of a norm squared in Eq. (p^. For t = one must have Wi(0) = /(O) 
since for short times any photon must be the first. This identity is easily checked 
by means of Eqs. (0), (0) and 

For equal dipole moments and without laser the conditional Hamiltonian is 
diagonal in the Dicke basis. A ± ReC describes the decay rates of \s) and \a) 
to \g), while ±h ImC can be viewed as a level shift. The state |e) can decay 
to both \s) and \a), with respective decay rates A ± ReC. This also follows 
from the Bloch equations and is indicated in Fig. |1|. From this the well-known 
fact follows that two atoms with dipole interaction can decay faster or slower 
than two independent atoms (super- and sub-radiance [^). When r — 0, ReC 
approaches A so that \a) can no longer decay while |s) decays with 2 A. 

Trajectories and Bloch equations. Starting at t = with a pure state, the state 
develops according to f/cond until the first emission at some time ti, determined 
from wi in Eq. (P^l). Then the state is reset according to Eq. ( [T7| ) to a new 
density matrix (which has to be normalized), and so on. 

The decomposition of 7^(p) in Eq. (^) allows one, however, to work solely 
with pure states which is numerically much more efficient. One can start with 
a pure state \ip), develop it with t/cond until ti to the (non-normalized) \i/j(ti)), 
reset to one of the pure states R±\ip{ti)) /\\ ■ \\ with relative probabilities given by 
the factors A ± ||Ci2 + C21I appearing in Eq. (pOD, and so on. The waiting time 
distributions are not changed by this procedure. 

Quite generally the ensemble of such trajectories yields the Bloch equations 
[ P7[] . With the reset matrix this is easily seen as follows. If an ensemble of 
systems of two two-level atoms has a density matrix p{t) at time t then at time 
t + At one has two sub-ensembles, one with a photon emission, the other with 
none. The former has relative size tr 7?.(p(t))At, by the remark after Eq. ([171) , 
while the latter is obtained by means of f/cond (^ + At, t) = 1a— ^-ffcondAt. This 
immediately gives 

P=-^ WcondP - pHlj + 7^(p) . (24) 



13 



Inserting from Eqs. (|2T| ) and ( ]T7| ) one obtains the Bloch equations for two two- 
level atoms. They agree with those derived by Agarwal 0. From this expression 
it is evident that -ffcond or the reset matrix can be immediately determined if the 
Bloch equations and the reset matrix or, respectively, -ffcond are explicitly known. 
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